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ABSTRACT 

Conventional  reliability-based  design  optimization  (RBDO) 
uses  the  means  of  input  random  variables  as  its  design  variables; 
and  the  standard  deviations  (STDEVs)  of  the  random  variables 
are  fixed  constants.  However,  the  fixed  STDEVs  may  not 
correctly  represent  certain  RBDO  problems  well,  especially 
when  a  specified  tolerance  of  the  input  random  variable  is 
presented  as  a  percentage  of  the  mean  value.  For  this  kind  of 
design  problem,  the  coefficients  of  variations  (CO Vs)  of  the 
input  random  variables  should  be  fixed,  which  means  STDEVs 
are  not  fixed.  In  this  paper,  a  method  to  calculate  the  design 
sensitivity  of  probability  of  failure  for  RBDO  with  fixed  COV  is 
developed.  For  sampling-based  RBDO,  which  uses  Monte  Carlo 
simulation  for  reliability  analysis,  the  design  sensitivity  of  the 
probability  of  failure  is  derived  using  a  first-order  score  function. 
The  score  function  contains  the  effect  of  the  change  in  the 
STDEV  in  addition  to  the  change  in  the  mean.  As  copulas  are 
used  for  the  design  sensitivity,  correlated  input  random  variables 
also  can  be  used  for  RBDO  with  fixed  COV.  Moreover,  the 
design  sensitivity  can  be  calculated  efficiently  during  the 
evaluation  of  the  probability  of  failure.  Using  a  mathematical 
example,  the  accuracy  and  efficiency  of  the  developed  method 
are  verified.  The  RBDO  result  for  mathematical  and  physical 
problems  indicates  that  the  developed  method  provides  accurate 
design  sensitivity  in  the  optimization  process. 
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NOMENCLATURE 

RBDO 

reliability-based  design  optimization 

STDEV 

standard  deviation 

COV 

coefficient  of  variation 

PDF 

probability  density  function 

CDF 

cumulative  distribution  function 

MCS 

Monte  Carlo  simulation 

FDM 

finite  difference  method 

N 

number  of  random  variables 

NC 

number  of  constraints 

NDV 

number  of  design  variables 

xt,  X 

random  variable,  random  variable  vector,  X  = 

[*1 . XNf 

Xu  X 

realization  of  Xt  and  X,  x  =  [x1, ... ,  xN]T 

cf  j,  d 

design  variable,  design  variable  vector,  d  = 

[d1, ... ,  dNDV\ 

Gj  (X) 

/-tli  performance  measure,  Gj  (X)  >  0  indicates 

failure. 

/x(x) 

joint  PDF  of  X 

fxSxi.dobd 

marginal  PDF  of  random  variable  XL 

FXl(.xi>  ai,bi) 

marginal  CDF  of  random  variable  Xt 

lii.Oi.COVi 

mean,  STDEV  and  COV  of  random  variable  Xt 

at, 

two  parameters  for  marginal  PDF  and  CDF  of 
random  variable  Xt 

(first-order)  score  function  of  design  variable. 

c(l -  91n/x(x) 

■V  w  -  ddi 
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(first-order)  score 
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function  of  STDEV, 


1.  INTRODUCTION 

Reliability-based  design  optimization  (RBDO)  has  been 
developed  using  the  mean  of  the  input  random  variable  as  the 
design  variable.  At  the  same  time,  the  standard  deviation 
(STDEV)  of  the  input  random  variable  is  invariant  in  the 
conventional  RBDO  process.  Many  problems  have  been  studied 
and  solved  assuming  the  STDEV  is  fixed  in  the  RBDO  process 
[1-10].  However,  the  fixed  STDEV  of  an  input  random  variable 
may  not  represent  physical  design  problems  appropriately.  In  the 
physical  design  problems,  the  tolerance  of  input  random 
variables  is  given.  The  tolerance  indicates  the  range  within 
which  the  input  random  variable  can  vary,  and  it  is  usually  a  fixed 
percentage  of  the  mean  of  the  input  random  variable,  which  is 
the  design  variable  in  RBDO.  Therefore,  the  range  changes  as 
the  design  variable  changes  during  the  RBDO  process  according 
to  the  percentage.  The  range  can  be  interpreted  as  the  STDEV  of 
the  input  random  variable;  hence,  STDEV  should  be  considered 
to  change  during  the  RBDO  process.  In  this  case,  the  coefficient 
of  variation  (COV),  which  is  the  ratio  of  STDEV  to  mean,  should 
be  fixed  in  the  RBDO  process  to  represent  the  given  tolerance. 

Fixing  the  COV  implies  that  the  STDEV  of  the  input  random 
variable  is  linearly  changing  as  the  design  variable  changes;  and 
a  varying  STDEV  indicates  that  the  uncertainty  of  the  input 
random  variable  is  changing  as  design  iteration  proceeds.  If  the 
uncertainty  of  the  input  random  variable  changes  in  each  design 
iteration,  the  RBDO  optimum  design  could  become  a  moving 
target  problem  unless  accurate  design  sensitivity  of  probability 
of  failure,  which  is  the  probabilistic  constraints  in  RBDO,  is 
provided.  In  sensitivity-based  RBDO,  such  as  the  first-order 
reliability  method  [1,  2,  11,  12],  the  second-order  reliability 
method  [13,  14],  and  the  dimension  reduction  method  [3,  15],  the 
most  probable  point  (MPP)  needs  to  be  found.  However,  the 
accurate  position  of  MPP,  where  the  design  sensitivity  is 
calculated  [16],  cannot  be  easily  obtained  because  the  position 
of  the  MPP  would  change  significantly  compared  to  the  case 
when  the  STDEV  is  fixed.  The  reason  for  this  is  that  the 
transformation  fromX-space  (actual  design  variable  space)  to  U- 
space  (independent  standard  normal  space),  which  is  necessary 
to  find  the  MPP,  is  affected  by  both  the  mean  and  the  STDEV, 
which  is  being  changed.  Hence,  an  approximated  RBDO  process 
such  as  the  Sequential  Optimization  and  Reliability  Assessment 
(SORA)  method  [6],  which  does  not  find  an  accurate  MPP  point 
in  the  design  process,  would  suffer  difficulty.  That  is  why  a 
specific  method  was  developed  for  RBDO  with  fixed  COV  using 
SORA  [17].  The  sampling-based  method  does  not  use  MPP; 
therefore,  the  same  problem  would  not  be  encountered  in 
sampling-based  RBDO.  However,  the  design  sensitivity  method 
for  sampling-based  RBDO  with  fixed  COV  has  not  been 
developed  yet. 


In  this  paper,  the  design  sensitivity  of  probability  of  failure 
for  sampling-based  RBDO  with  fixed  COV  is  developed.  First, 
the  sampling-based  RBDO  method  is  briefly  revisited,  and  then 
the  score  functions  of  the  mean  and  the  STDEV  of  the  input 
random  variable  are  derived.  In  addition  to  marginal  PDF, 
copulas  [18,  19]  are  used  in  the  derivation  of  score  function  to 
consider  correlated  input  random  variables.  Next,  the  score 
functions  are  used  to  calculate  design  sensitivity  for  RBDO  with 
fixed  COV.  Finally,  the  proposed  design  sensitivity  is  verified 
using  numerical  examples  of  RBDO  with  fixed  COV. 

2.  REVIEW  OF  SAMPLING-BASED  RBDO 

A  general  RBDO  problem  is  expressed  as 


minimize  cost(d) 


subject  to 


P[Gj(X)  >  0]  <  P^r,  j  =  1, ...  ,NC 
dL  <  d  <  du,  d  £  Mndv,  and  X  £ 


(1) 


where  X  is  the  V-dimensional  random  variable  vector,  d  is  the 
MDF-dimensional  design  variable  vector,  dL  is  the  lower  design 
bound,  d,;  is  the  upper  design  bound,  P^r  is  the  target 
probability  of  failure  for  the  j- th  performance  measure  Gy(X), 
and  NC  is  the  number  of  constraints.  In  sampling-based  RBDO, 
the  probabilistic  constraint  -  probability  of  failure  -  is  estimated 
using  the  Monte  Carlo  simulation  (MCS)  method  as  [4,  5,  20] 


Pf,  =  P[GjQQ  >  0]  =  f  V(x)/x(x)  dx 
Jnx  > 

nMCS 


(2) 


<?=1 


where  nMCS  is  the  number  of  MCS  samples,  x(£?)  is  the  g-th 
realization  (MCS  sample)  of  X,  a  failure  set  f \F.  is  defined  as 
PlF.  =  (x:  Gj  (x)  >  0],  and  the  indicator  function  is  defined  as 


x  £  nF., 
otherwise. 


(3) 


In  Eq.  (2),  it  can  be  seen  that  sampling-based  RBDO  does 
not  require  the  sensitivity  (gradient)  of  the  performance  measure 
Gj  (X)  to  evaluate  the  probability  of  failure.  This  is  a  significant 
merit  for  engineering  problems  in  which  it  is  difficult  to  obtain 
the  gradient  (e.g.,  explosion  analysis)  because  RBDO  can  be 
accomplished  without  the  gradient  information  using  sampling- 
based  RBDO.  On  the  other  hand,  demerit  of  sampling-based 
RBDO  is  that  the  accurate  estimation  of  probability  of  failure 
depends  on  the  number  of  MCS  samples  nMCS.  The  existence 
of  many  nMCS  indicates  accurate  estimation;  however,  it  also 
requires  many  evaluations  of  the  performance  measure  G;(X), 
which  could  be  computationally  expensive.  Hence,  a  surrogate 
model  is  used  to  calculate  the  probability  of  failure  in  Eq.  (2)  [4, 
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5].  In  this  paper,  the  surrogate  model  of  performance  measure  is 
assumed  to  be  available. 

The  probability  of  failure  expression  in  Eq.  (2)  requires  the 
input  joint  PDF  of  fx (x) .  Also,  MCS  samples  to  calculate  the 
probability  of  failure  should  be  generated  according  to  the  input 
joint  PDF.  In  this  paper,  the  input  joint  PDF  is  defined  using  the 
marginal  PDF  of  each  input  random  variable  and  copula  for  each 
correlated  pair  of  input  random  variables.  Copula  can  connect 
two  or  more  marginal  PDFs  of  random  variables  to  represent 
statistical  correlation  between  the  random  variables  and  can 
generate  a  joint  PDF  using  the  copula  density  function.  In  this 
paper,  statistical  correlation  between  two  random  variables  is 
considered.  When  the  A-th  correlated  random  variable  pair  Aj¬ 
ar)  d  Xj  are  correlated  with  the  copula  density  function  ck,  the 
joint  PDF  ofYj-  and  Xj  can  be  represented  as  [18,  19,  23] 


variables  as  well  [22].  In  addition,  a  score  function  method  to 
calculate  the  design  sensitivity  of  the  probability  of  failure  has 
been  developed  for  sampling-based  RBDO  with  fixed  STDEV 
for  correlated  input  random  variables  [4,  5].  In  this  paper,  the 
design  sensitivity  using  the  score  function  is  further  expanded  to 
apply  to  sampling-based  RBDO  for  correlated  input  random 
variables  with  fixed  COV. 

The  design  sensitivity  of  the  probability  of  failure  can  be 
obtained  by  taking  the  derivative  of  the  probability  of  failure 
expression  in  Eq.  (2)  with  respect  to  a  design  variable  dt  as 


fxtXj{xt,X})  (4 

=  fxSXi.  a-i,  bt )  fXj{xj ;  ap  bj)  ck(u,  v ;  6k) 

where  fXl  (xt;  a b[)  is  the  marginal  PDF  of  the  / -th  input  random 
variable  Xh  at  and  bL  are  two  parameters  of  the  marginal 
distribution,  ck(u,v;  0k )  is  the  copula  density  function  of  the  A- 
th  correlated  pair,  u  and  v  are  the  cumulative  distribution 
function  (CDF)  values  of  the  correlated  pair,  and  9k  is  the 
correlation  coefficient  for  the  copula.  In  this  paper,  it  is  assumed 
that  each  random  variable  has  marginal  distribution  with  two 
parameters,  as  is  usually  assumed  in  the  RBDO  process.  Using 
the  relationship  in  Eq.  (4),  the  input  joint  PDF  /x(x)  can  be 
expressed  as  [18,  19,  23] 

N  M 

/x(x)  =  Y\  /*;(**;  ai’  ht>  n  CM- v'  dj)  (5) 

i= 1  7=1 

where  N  is  the  number  of  random  variables,  and  M  is  the  number 
of  correlated  pairs.  It  is  noted  that  the  mean  and  the  STDEV 
ak  of  the  A-th  random  variable  Xk  are  only  related  to  the  marginal 
PDF  fXk(xk,  afc,bfc)  and  copula  density  function  cm(u,  v;  6m) 
when  Xk  is  in  the  m- th  correlated  pair.  The  input  joint  PDF 
defined  in  Eq.  (5)  can  be  used  to  calculate  the  probability  of 
failure  in  Eq.  (2)  and  to  obtain  the  design  sensitivity  of  the 
probability  of  failure  in  the  following  section. 

3.  DESIGN  SENSITIVITY  AND  SCORE  FUNCTIONS 

For  an  effective  optimization  process,  an  accurate  and 
efficient  design  sensitivity  needs  to  be  provided,  especially  for 
RBDO  with  fixed  COV.  The  design  sensitivity  of  the  probability 
of  failure  for  sampling-based  RBDO  can  be  obtained  using  the 
score  function.  The  score  function  of  a  variable  is  a  partial 
derivative  of  the  log  function  (log-derivative)  of  the  input  joint 
PDF  in  Eq.  (5)  with  respect  to  the  variable.  The  score  function 
to  estimate  the  stochastic  sensitivity  of  the  general  functions  of 
independent  random  variables  was  developed  first  [21].  The 
method  was  applied  to  functions  of  correlated  Gaussian  random 


In  Eq.  (6),  it  can  be  seen  that  a  new  term,  the  log-derivative  of 
the  input  joint  PDF  /x(x)  with  respect  to  the  design  variable  dh 
is  present  in  addition  to  the  probability  of  failure  expression  in 
Eq.  (2).  This  additional  term  is  the  first-order  score  function  of 
the  design  variable  dt  as 


d  ln/x(x) 
ddt 


(7) 


It  is  noted  that  the  first-order  score  function  in  Eq.  (7)  does  not 
include  the  sensitivity  (gradient)  of  the  performance  measure 
Gj  (X) .  Therefore,  the  design  sensitivity  of  the  probability  of 
failure  in  Eq.  (6)  can  be  calculated  without  the  gradient  of  G;-(X). 
Using  the  MCS  method  and  the  first-order  score  function,  the 
design  sensitivity  in  Eq.  (6)  can  be  estimated  similar  to  Eq.  (2) 
as  [4,  5,  20] 

PFj  =  [  w . (x)  (x)  fx (x)  dx 

Jpix  1 

nMCS  (8) 

-ass  I  Sl-^1 

?= 1 


Let  the  design  variable  dt  be  /rfc,  which  is  the  mean  of  the  A- 
th  random  variable  Xk.  Also,  let  ak  and  COVk  be  the  STDEV  and 
COV  of  the  random  variable  Xk.  Then  the  first-order  score 
function  in  Eq.  (7)  can  be  further  derived  as 

^f(x) 

_p“(x)  if  Gk  is  fixed,  (9) 

(s^Cx)  +  COVk  5®(x)  if:  COVk  is  fixed. 

where  sj£  (x)  =  dln/x(-x>  and  5®  (x)  =  dln/x(-x)  are  first-order 

it  dak 

score  functions  of  mean  and  STDEV,  respectively.  Hence,  to 
evaluate  the  design  sensitivity  in  Eq.  (6),  the  score-functions 
S^k  (x)  and  S„k  (x)  are  required.  As  explained  earlier,  /rk  and 
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<jk  are  related  only  to  marginal  PDF  fXk(xk',  ak<  bk)  and  copula 

density  function  cm  (u,  v, ;  0)n)  in  the  joint  PDF  /x(x)  in  Eq.  (5). 

Cl)  Cl) 

Therefore,  the  score  functions  (x)  and  Sak (x)  can  be 
obtained  as 


and 


d  ln/x(x) 

3/tfc 

d\nfXk{xk-,ak,bk )  |  d  In  cm(u,  v,  6>;n) 
9Pfc  Sylifc 


(10) 


dln/x(x) 

dak 

d\nfXk{xk-,ak,bk)  |  d  In  cm(u,  v\  9?n) 
dak  dak 


(11) 


In  Eqs.  (8)  and  (9),  it  can  be  seen  that  log-derivatives  of 
fxk(.xk'’  ak>  bk)  and  cm (M’  v’  J  ^?n)  with  respect  to  fik  and  ak  have 
to  be  derived.  Flere,  it  is  noted  that  the  log-derivative  of  the 
copula  density  function  can  be  ignored  when  the  input  random 
variable  Xk  is  statistically  independent. 

3.1  Log-derivatives  of  marginal  PDF  with  respect  to  mean  and 
STDEV 

Seven  widely  used  marginal  distribution  types  are 
considered  in  this  paper.  Their  PDFs  and  CDFs  are  shown  in 
Appendix,  Table  A.l.  It  can  be  seen  that  they  have  two 


parameters  a  and  b  in  their  analytical  forms.  The  two  parameters 
can  uniquely  determine  the  mean  and  the  STDEV  of  the  input 
random  variable  and  vice  versa,  as  shown  in  Appendix,  Table 
A. 2.  Therefore,  the  log-derivatives  of  the  marginal  PDF  in  Eqs. 
(10)  and  (11)  can  be  represented  using  the  two  parameters  a  and 
b  as 


and 


d\nfXk{xk,ak,bk) 

dln/xfc(*fc;%A) 

dak 


d\nfXk  dak  |  d\nfXk  dbk 
dak  dfik  dbk  d/rk 


(12) 


d\nfXk  dak  |  d\nfXk  dbk 
dak  dak  dbk  dak 


(13) 


where  ak  and  bk  are  the  two  parameters  of  the  input  random 
variable  Xk.  The  right  sides  of  Eqs.  (12)  and  (13)  are  composed 
of  the  log-derivatives  of  the  marginal  PDF  with  respect  to  the 
two  parameters  and  their  derivatives  to  the  mean  and  STDEV. 
Using  the  analytical  expression  of  the  marginal  PDF  in  Table 
A.l,  the  log-derivatives  of  the  marginal  PDF  with  respect  to  the 
two  parameters  a  and  b  are  obtained  as  shown  in  Table  1.  Using 
the  relationship  in  Table  A.2,  derivatives  of  a  and  b  with  respect 
to  the  mean  q  and  the  STDEV  a  are  also  obtained  as  shown  in 
Table  2.  By  substituting  the  derivative  terms  in  Tables  1  and  2 
for  Eqs.  (12)  and  (13),  the  log-derivatives  of  the  marginal  PDF 
with  respect  to  the  mean  and  STDEV  in  Eqs.  (10)  and  (11)  can 
be  obtained. 


Table  1  Log-derivatives  of  Marginal  PDF  with  Respect  to  Parameters 


Distribution 

type 

5  In  fXi 

ddi 

5  In  fXi 
dbi 

Normal 

Xi  Cl  i 

bf 

1  (Xi  -  Ui)2 

bi  bf 

Lognormal 

In  xt  —  aL 

bf 

1  (lnXj  —  af)2 
bi  +  bf 

Weibull 

M-i+pfl 

CL[  'Q-i' 

1  (xi 

—  +  1  —  —  (In  Xi  -  In  af) 

bi  V&j/ 

Gumbel 


Gamma 


Extreme 


Extreme  Type 

II 


1 


lnXj  —  ^(a;)  —  In  bi* 


1 

— I-  (In  bt  —  In  Xj) 

G-i 


l 

bi 


bf 


1 

bi 


!/)(•):  Digamma  function,  i p(s)  =  r'(s)/F(s) 
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Table  2  Derivatives  of  Parameters  with  Respect  to  Mean  and  STDEV 


Distribution 

type 

da.i 

d\ii 

dbi 

dUi 

da.i 

doi 

dbi 

d<Ji 

Normal 

1 

0 

0 

1 

Lognormal 

1  (l  1  ^ 

of 

A 

A 

A  V  A2  +  of) 

bi  Mi  (A2  +  of) 

A2  +  of 

bi  (A2  +  of) 

Weibull 

'dai/d/ii 

dbi/d^ 

=  A"1 

1 

.A 

dai/dui 

dbi/diJi 

=  A"1 

0 

A 

where  A  = 

r(i  + 1  /bi)  -  ±  r(i  + 1 /bi)  1/;(1  +  l /bi) 

°i 

a,  r(l  +  2/bi)  -  S  r(i  +  2/bi)  <K  1  +  2/bi) 

°i 

Gumbel 

1 

0 

ai 

-3.4632—1- 

7 Tzbi 

6  <7j 
n2bi 

Gamma 

2 

bi 

1 

CLi 

2-,/a 

bi 

2 

■yj  Cl  i 

Extreme 

1 

0 

0.5772-06 

V6 

n 

n 

Extreme 
Type  11 

dai/d/ii 

dbi/dfii. 

=  B1 

l 

A. 

dai/dui 

dbi/dai_ 

=  B-1 

O' 

.A. 

where  B  = 

'5  -  Va)  0(i  -  1/A)  r(i  -  1/ a)  ’ 

ai 

%  r(l  -  2 /CLi)  0(1  -  2/a )  bi  r(l  -  2 /a/) 

_i 

3.2  Log-derivatives  of  copula  density  function  with  respect  to 
mean  and  STDEV 

The  second  terms  in  the  right  sides  ofEqs.  (10)  and  (11)  are 
the  log-derivatives  of  the  copula  density  function.  As  explained 
earlier,  the  log-derivatives  of  the  copula  density  function  can  be 
ignored  when  Xk  is  an  independent  random  variable  because  the 
copula  density  function  for  independent  random  variables  is  one 
and  its  log-derivative  is  zero.  If  is  statistically  correlated  and 
u  =  Fxk  ( xk,ak,bk )  is  the  CDF  value  of.ri-,  the  log-derivatives 
of  the  copula  density  function  become 


d  In  cm  ( u ,  v ;  6m) 

dlncm 

du 

(14) 

d^k 

du 

d\ik 

d\ncm(u,v;6m) 

d  In  cm 

du 

(15) 

dak 

du 

dak' 

For  various  copula  types,  the  log-derivatives  of  the  copula 
density  function  with  respect  to  u  (dlncm  /du),  in  the  right 
sides  ofEqs.  (14)  and  (15),  are  well  summarized  in  Ref.  4.  The 
other  terms  du  /d\ik  and  du  /dok  are  obtained  in  this  paper.  In 
Table  A.  1,  it  can  be  seen  that  the  marginal  CDF  of  the  input 
random  variable  Xk,u  —  Fxk  ( xk,  ak,  bk)  is  also  a  function  of  the 


two  parameters  ak  and  bk.  Therefore,  the  derivatives  of  CDF  u 
can  be  expressed  using  the  two  parameters  as 


and 


du  du  dak  du  dbk 
d/ik  dak  dfik  dbk  djik 


du  du  dak  du  dbk 
dak  dak  dak  dbk  dak' 


(16) 

(17) 


In  Eqs.  (16)  and  (17),  the  derivative  of  u  with  respect  to  the 
two  parameters  ak  and  bk  are  necessary.  Using  the  analytical 
representations  of  the  marginal  CDFs  in  Table  A.l,  the 
derivatives  of  the  seven  marginal  CDFs  are  derived  and 
summarized  in  Table  3.  Using  Ref.  4  and  Table  3,  the  log- 
derivatives  of  the  copula  density  function  in  Eqs.  (14)  and  (15) 
are  available.  Then,  all  terms  for  Eq.  (9)  are  obtained,  and  the 
design  sensitivity  of  the  probability  of  failure  for  sampling-based 
RBDO  with  fixed  COV  can  be  calculated  using  Eq.  (8). 


3.3  Score-functions  for  sampling-based  RBDO  with  fixed  COY 
Using  the  derivatives  obtained  in  previous  sections,  the 
design  sensitivity  of  probability  of  failure  for  sampling-based 
RBDO  with  fixed  COV  can  be  obtained.  Unlike  the  case  in 
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Table  3  Derivative  of  CDF  with  Respect  to  Parameters 


Distribution 

type 

du 

dat 

du 

~dbt 

Normal 

1  (Xi-  a  a 

bi  A  6;  J 

-  tf  n  ^  ) 

Lognormal 

1  /In  Xi  —  aA 
b^{  ^  ) 

lnXj  —  a.i  /In  Xi  -  ao 
b?  0  {  ^  ) 

Weibull 

bi  /Xi\bi 

(1  -«) 

CLi  \CLi/ 

(Xt \bi 

(In  xt  —  In  a,  )  ^ — J  (1  —  u) 

Gumbel 

1  /  Xi-CLA 

biexp{  bi  )u 

Xi  —  CLi  (  xi  —  ai\ 

bf  expl  ^  )u 

Gamma 

r(ai)xp(ai')  —  \nXi  T 

—  Xj7’(3,aj,Xj/h,)t 

cT 

S' 

i 

tu 

1 

Extreme 

1  (Xi  -  a,\ 

i,,exp(  b,  )(“ 

Xi  —  a,-  (Xi  —  CLj\ 

bf  exp(  ^  j*”-1) 

Extreme  Type 

II 

(bi\ai 

—  (In  bt  —  In  x^  J  u 

bt  \XiJ 

which  STDEV  is  fixed,  the  design  sensitivity  for  RBDO  with 
fixed  COV  requires  all  terms  in  Eqs.  (10)  to  (17)  for  its 
calculation.  Therefore,  the  design  sensitivity  expression  for 
RBDO  with  fixed  COV  would  be  more  involved.  However,  for 
some  marginal  distribution  types,  the  expression  could  be 
simpler. 

The  second  expression  on  the  right  side  of  Eq.  (9)  is  the  first- 
order  score  function  of  design  variable  d,  for  RBDO  with  fixed 
COV.  When  Eqs.  (10)  to  (17)  are  substituted  for  the  expression, 
the  score  function  becomes 

SdfOO 

/dhi  fXk  d  In  cm  du\/dak  dak\ 

y  dak  du  dak)  \d(ik  k  dak)  (18) 
/din  fXk  d  In  cm  du\/dbk  dbk\ 

y  dbk  du  dbk)\d(ik  k  dak) 


= 

din fXk  d In cm  du  \  / dak  ^  dajA  (20) 
dak  du  dakJ\d(ik  k  dak) 


Also,  for  the  Weibull  distribution,  Eq.  (19)  holds; 
furthermore,  the  relationship  in  Eq.  (21),  which  can  be  derived 
using  the  derivatives  in  Table  2,  holds  as  well. 


dak  dak  1 

_ _  _i_  rnv  _ _  — _ 

dMfc  k  dfffc  f(l  +  l/bk) 


(21) 


Hence,  the  score  function  in  Eq.  (18)  can  be  expressed  in  a 
simpler  form  for  the  Weibull  distribution  as 


/d  In  fXk  d  In  cm  du  \  1 

y  dak  +  du  dak)  f(l  +  l/bk)' 


(22) 


In  Table  2,  the  derivatives  of  the  two  parameters  a  and  b  with 
respect  to  mean  and  STDEV  have  been  shown.  It  can  be  found 
that  the  following  relationship 


dbk  dbk 

-±+COVk-±=  0 
dMfc  d<rk 


(19) 


holds  for  the  Lognormal  distribution.  Therefore,  Eq.  (18)  for  the 
Lognormal  distribution  becomes  a  simpler  form  as 


The  Extreme  type-II  distribution  also  has  special 
relationships  when  COV  is  fixed  as 


and 


dak 

d\ik 


+  COVk 


dak 

dak 


■=  0 


dbk  dbk  1 

■JT±  +  COVk-±  =  — - — -. 

d/rk  dak  f(l  -  l/ak) 


(23) 

(24) 


t  r(.,.):  Upper  incomplete  gamma  function,  r(s,x)  =  f  ts  1e  1  dt , 
T  A  case  of  Meijer  G-function.  See  Ref.  24. 
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Finally,  the  Extreme  type-II  distribution  has  a  simpler  form  of 
Eq.  (18)  as 


/din  fXk  d  In  cm  du  \  1 

y  dbk  +  du  dbk)  T(1  —  l/ak)' 


(25) 


Because  the  score  function  S^.  (x)  is  derived  analytically, 
the  design  sensitivity  of  probability  of  failure  in  Eq.  (8)  can  be 
calculated  along  with  calculation  of  the  probability  of  failure  in 
Eq.  (2)  with  small  additional  effort.  Hence,  the  design  sensitivity 
for  RBDO  with  fixed  COV  can  be  evaluated  efficiently  using  the 
score  functions  during  the  calculation  of  the  probability  of 
failure. 


4.  NUMERICAL  EXAMPLES 

The  effectiveness  and  efficiency  of  the  developed  method 
for  the  design  sensitivity  of  the  probability  of  failure  are  tested 
using  numerical  examples  in  this  section.  The  developed  design 
sensitivity  is  compared  to  the  sensitivity  using  the  finite 
difference  method  (FDM)  to  verify  its  accuracy.  The  FDM 
results  are  obtained  after  a  number  of  trial-and  error  process  in 
the  determination  of  appropriate  step  size  to  yield  accurate 
design  sensitivity.  Then,  using  the  developed  design  sensitivity, 
RBDO  with  fixed  COV  is  performed  for  a  2-D  mathematical 
problem  and  a  1 4-D  engineering  problem  to  understand  ho  w  well 
the  developed  design  sensitivity  method  works  for  the  RBDO 
process. 


4.1  Accuracy  check 

The  accuracy  of  the  developed  design  sensitivity  of 
probability  of  failure  is  verified  using  a  highly  nonlinear 
performance  measure 

G(X)  =  0.7361  +  (F  —  6)2  +  (F  —  6)3  . 

-0.6  X  (F  -  6)4  +  Z  (  t>) 

where 

m  _  ro.9063  0.4226  i 

Yz\  ~  1.0.4226  -0.9063.1  U2J 

and  G(X)  >  0  indicates  failure.  For  thorough  testing,  five 
different  input  models,  which  include  all  marginal  distribution 
and  copula  types  dealt  with  in  the  previous  sections  are  used  for 
the  testing.  The  design  variables  d ±  and  d2  are  the  means  of  the 
input  random  variables  X1  and  X2,  respectively.  In  the  design 
sensitivity  calculation,  the  CO  Vs  of  both  random  variables  are 
fixed. 

The  accuracy  of  the  developed  design  sensitivity  of  the 
probability  of  failure  is  verified  by  comparing  it  with  the  forward 
FDM  design  sensitivity.  The  result  of  the  developed  design 
sensitivity  and  the  FDM  design  sensitivity  is  summarized  in 
Table  5.  The  ratio  of  the  developed  sensitivity  to  the  FDM 
sensitivity  is  listed  in  the  ‘Accuracy’  column  of  the  table.  The 
ratio  varies  from  97.1%  to  102.6%,  which  indicates  that  the 
developed  sensitivity  agrees  with  the  FDM  design  sensitivity. 
Therefore,  it  can  be  concluded  that  the  developed  design 
sensitivity  is  accurate. 


Table  4  Input  Models  for  Accuracy  Test 


V, 

x2 

Copula 

Type  Mean  STDEV  COV 

Type  Mean  STDEV  COV 

Type  Tau 

A 

Weibull 

Normal 

Clayton  0.3 

B 

Gumbel 

Lognormal 

Frank  0.3 

C 

Gamma  4  0.5  0.125 

Extreme  2.5  0.5  0.2 

FGM  0.2 

D 

Extreme 

Normal 

Gaussian  0.3 

E 

Extreme  II 

Lognormal 

AMH  0.33 

Table  5  Comparison  of  FDM  and  Developed  Design  Sensitivity 


Input 

d  i 

di 

model 

Perturb 

FDM 

Developed 

Accuracy 

Perturb 

FDM 

Developed 

Accuracy 

A 

0.1% 

0.600 

0.597 

99.5% 

0.1% 

-0.135 

-0.133 

98.5% 

B 

0.1% 

0.346 

0.338 

97.7% 

1.0% 

-0.107 

-0.106 

99.1% 

C 

0.1% 

0.458 

0.469 

102.4% 

1.0% 

-0.174 

-0.169 

97.1% 

D 

0.1% 

0.639 

0.638 

99.8% 

0.5% 

-0.144 

-0.147 

102.1% 

E 

0.2% 

0.309 

0.306 

99.0% 

0.5% 

-0.077 

-0.079 

102.6% 
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In  this  example,  the  FDM  design  sensitivity  of  the 
probability  of  failure  is  used  as  a  reference  value  to  check  the 
accuracy  of  the  developed  design  sensitivity.  To  obtain  the 
accurate  FDM  design  sensitivity,  two  conditions  -  appropriate 
perturbation  size  and  highly  accurate  probability  of  failure  at 
both  perturbed  and  original  design  -  should  be  satisfied. 
However,  obtaining  appropriate  perturbation  size  is  not  an  easy 
task.  A  large  perturbation  size  cannot  provide  accurate  design 
sensitivity  of  a  nonlinear  function,  while  a  small  perturbation 
size  can  suffer  numerical  noise.  In  this  example,  perturbation 
sizes  of  0.1%,  0.2%,  0.5%  and  1%  of  the  design  variable  value 
are  used,  as  shown  in  Table  5.  For  different  cases  and  different 
design  variables,  perturbation  size  should  be  carefully  chosen; 
hence,  several  perturbation  sizes  are  tested,  and  the  perturbation 
size  that  shows  a  stable  design  sensitivity  value  is  selected.  An 
appropriate  perturbation  size  at  one  design  configuration  may 
not  be  applicable  at  the  other  design  configuration  because  the 
nonlinearity  of  the  probability  of  failure  changes  as  the  design 
changes.  Hence,  it  is  not  easy  to  use  the  FDM  design  sensitivity 
in  the  optimization  process.  While  the  developed  design 
sensitivity  method  does  not  depend  on  the  perturbation  size,  it 
shows  the  same  accuracy  as  the  FDM  design  sensitivity;  hence, 
it  is  more  suitable  for  the  optimization. 

The  second  condition  -  accurate  probability  of  failure  -  is 
also  important  in  order  to  obtain  the  correct  design  sensitivity  of 
the  probability  of  failure  using  the  FDM.  In  Table  5,  the  FDM 
design  sensitivity  to  d1  in  input  model  A  is  0.600.  To  obtain  this 
value,  the  finite  difference  0.00240  between  two  probabilities  of 
failure,  0.13726  at  the  perturbed  design  and  0.13486  at  the 
original  design,  is  divided  by  the  perturbation  size  0.004,  which 
is  0.1%  of  d±.  It  is  noted  that  the  significant  digit  in  the  two 
probabilities  of  failure  should  be  larger  than  five  to  obtain  three 
significant  digits  in  the  finite  difference  value.  If  an  estimated 
value  of  probability  of  failure  is  needed,  three  significant  digits 
might  be  enough.  However,  because  the  finite  difference  value 
is  needed  to  calculate  the  FDM  design  sensitivity,  more  accurate 
probabilities  of  failure  are  required  to  have  enough  significant 
digits  in  the  finite  difference  value  calculation.  For  this  reason, 
in  this  example,  200  million  MCS  samples  were  required  to 
accurately  calculate  the  probabilities  of  failure  at  the  original  and 
perturbed  design;  hence,  a  total  of  400  million  samples  are  used 
to  calculate  one  value  of  FDM  design  sensitivity.  On  the  other 
hand,  the  developed  method  evaluates  the  design  sensitivity 
without  calculating  the  finite  difference  value.  Hence,  we  can 
use  much  smaller  number  of  MCS  samples  to  calculate  an 
accurate  design  sensitivity.  In  this  example,  one  million  MCS 
samples,  which  are  0.25%  of  the  samples  for  the  FDM  design 
sensitivity,  are  used.  Hence,  it  can  be  seen  that  the  developed 
design  sensitivity  is  efficient  and  suitable  for  the  optimization 
process. 

4.2  2-D  mathematical  example 

To  find  how  well  the  developed  design  sensitivity  method 
works  in  the  optimization  process,  RBDO  is  performed  for  a  2- 
D  mathematical  example.  The  RBDO  formulation  is  shown  as 


.  .  .  „ ,  j,  (d1+d2-10)2  (d1-d2+10')2 

minimize  C(d)  = - 

v  '  30  120 

subject  to  P[G,(X)  >  0]  <  pj“r  =  2.275%,  (27) 

i  =  1,  ...,3,  dL  <  d  <  df 
d  E  R2  and  X  E  R2 

where  X=  [Xi,Z2]r,  d  =  [du  d2]r,  &L  =  [0,  Of,  du  =  [10,  10f, 
and  the  three  constraints  are  given  by 

X\X2 

1 - , 

20 

-1  +  (0.9063^  +  0.4226X2  -  6)2 
+  (0.9063X1  +  0.4226XZ  -  6)3 
—0.6(0. 90632^  +  0.4226XZ  -  6)4  1 

—  (— 0.4226Xl  +  0.9063X2), 

80 

1  - • 

X\  +  QX2  +  5 

The  contours  of  the  cost  function  C(d)  in  Eq.  (27)  and  limit 
states  of  the  three  constraints  in  Eq.  (28)  are  shown  in  Fig.  1. 


Gi(x)  = 
G2(X)  = 


G3(X)  = 


Figure  1  Cost  Function  and  Limit  State  Contours  for 
2-D  Mathematical  Example 

In  this  example,  two  different  input  models  are  used  for 
RBDO  as  shown  in  Table  6.  The  two  models  have  the  same 
means,  STDEVs,  COVs,  copula  type,  and  Kendall’s  tau,  while 
the  marginal  distribution  types  of  20  and  X2  are  different.  The 
two  design  variables  d\  and  ch  of  the  example  are  the  means  of 
X\  and  X2,  respectively,  and  the  COVs  of  A i  and  X2  are  fixed 
during  the  RBDO.  To  calculate  the  probabilities  of  failure  and 
the  design  sensitivities  for  the  probabilities  of  failure,  one 
million  MCS  samples  are  used  in  this  example. 
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Table  6  Input  Models  for  2-D  Mathematical  Example 


Input 

A, 

x2 

Copula 

model 

Type 

Mean 

STDEV 

COV 

Type 

Mean 

STDEV 

COV 

Type  Tau 

A 

Lognormal 

5.19 

0.7785 

0.15 

Normal 

0.74 

0.185 

0.25 

Clayton  0.5 

B 

Normal 

Gumbel 

Table  7  RBDO  Result  of  2-D  Mathematical  Example 


Input 

Optimum  design 

Iter. 

Reli. 

Cost 

Probability  of  failure  (%) 

model 

d\ 

di 

eri 

er2 

Anal. 

G\ 

G2 

g3 

A 

4.6650 

3.3192 

0.6998 

0.8298 

8 

9 

-1.2082 

2.2864 

2.2986 

1.5631 

B 

4.7610 

2.8072 

0.7141 

0.7018 

9 

10 

-1.3879 

2.2678 

2.2955 

0.7259 

(a)  Input  mode /  A  (b)  Input  mode!  B 


Figure  2  2-crContour  of  Input  Model  at  Initial  Design  and  Optimum  Design 


The  result  of  RBDO  using  the  developed  design  sensitivity 
for  RBDO  with  fixed  COV  is  shown  in  Table  7.  It  can  be  seen 
that  the  STDEVs  are  changed  as  the  design  variables  change. 
However,  the  CO  Vs  are  kept  during  the  optimization  process  to 
0.15  for  X\  and  0.25  for  A2.  Input  model  A  requires  eight  design 
iterations  and  nine  reliability  analyses.  This  means  that  only  one 
additional  line  search  is  needed  in  the  optimization  process,  as 
an  accurate  design  sensitivity  is  provided.  Similarly,  RBDO  of 
input  model  B  converges  in  nine  design  iterations  with  one 
additional  line  search,  which  indicates  the  effectiveness  of  the 
design  sensitivity  as  well.  In  both  cases,  a  2.275%  target 
probability  of  failure  is  satisfied  under  1%  relative  convergence 
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criteria  for  the  probabilistic  constraints.  Therefore,  both 
optimum  designs  satisfy  the  probabilistic  constraints. 

In  Fig.  2,  2-er  contours  of  the  input  model  at  initial  design  d° 
and  RBDO  optimum  design  dopt  are  shown.  As  the  limit  state 
functions  are  not  linear  and  some  input  random  variables  do  not 
follow  normal  distribution,  the  2-er  contours  at  RBDO  optimum 
designs  dopt  do  not  meet  the  limit  states  tangentially.  Note  that 
the  area  of  the  2-er  contours  are  dramatically  changed  from  the 
initial  design  to  the  optimum  design  because  the  STDEV  of  20 
changes  significantly  as  the  second  design  variable  di  changes. 
Especially  for  input  model  A,  the  2-er  contour  of  the  optimum 
design  fits  in  the  feasible  region  as  shown  in  Fig.  2(a),  so  a  slight 
movement  of  the  design  could  violate  one  or  more  probabilistic 
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constraints.  Hence,  unless  accurate  design  sensitivity  is 
provided,  optimum  design  will  not  be  easy  to  obtain  readily. 

The  same  RBDO  with  fixed  COV  has  been  carried  out  once 
again  using  inaccurate  design  sensitivity  by  ignoring  the  second 
term  COVk  (x)  in  the  second  design  sensitivity  expression  in 
Eq.  (9).  In  other  words,  the  effect  of  the  change  of  STDEV  is 
ignored.  The  RBDO  result  is  shown  in  Table  8.  In  both  cases, 
the  optimization  process  does  not  converge  to  the  optimum 
design  as  2.275%  target  probability  failure  cannot  be  achieved. 
In  addition,  the  number  of  reliability  analyses  is  over  18,  which 
is  almost  twice  that  of  previous  cases  when  an  accurate  design 
sensitivity  of  probability  of  failure  was  provided.  Therefore,  the 
STDEV  effect  in  the  design  sensitivity  calculation  is  necessary 
for  appropriate  and  efficient  convergence  to  the  correct  optimum 
design  of  RBDO  with  fixed  COV. 

4.3  14-D  engineering  example 

An  RBDO  has  been  performed  for  a  14-D  engineering 
problem.  This  example  is  a  car  noise,  vibration,  and  harshness 
(NVH)  and  crash-safety  problem.  There  are  11  constraints, 
which  include  full  frontal  impact,  40%  offset  frontal  impact,  and 
NVH,  as  shown  in  Table  9.  Surrogate  models  for  the  constraints 
and  cost  function  for  RBDO  has  been  provided  by  Ford  Motor 
Company.  Originally,  this  example  had  44  input  random 


variables,  which  means  it  was  a  44-D  RBDO  problem.  However, 
the  number  of  random  variables  has  been  reduced  to  14 
according  to  a  variable  screening  procedure  [25]. 

Input  random  variables  are  the  thickness  of  the  plates  in  a 
car  body.  They  have  Normal,  Lognormal,  Weibull,  and  Gamma 
distributions  as  marginal  distribution  types.  The  design  variables 
are  the  means  of  the  input  random  variables,  and  the  CO  Vs  are 
fixed  to  be  constants.  There  are  two  correlated  pairs  of  Xu  -  Xi 
and  X25  -  Xiu-  whose  correlations  are  represented  using  the 
Clayton  and  Gaussian  copulas,  respectively.  The  detailed 
information  of  the  input  model  is  shown  in  Table  10;  ds  indicates 
baseline  design,  which  is  the  initial  design  of  the  optimization. 

The  RBDO  formulation  of  the  14-D  engineering  example  is 
shown  in  Eq.  (28).  Cost  function  is  the  weight  of  the  car  body 
and  a  function  of  design  variable  vector  d,  not  input  random 
variable  vector  X.  Hence,  the  cost  function  is  a  deterministic 
function,  while  the  constraint  functions  are  a  random  function  of 
the  vector  X. 

minimize  C(d)  =  Weight(d) 

subject  to  P[G,(X)  >  Baseline J  <  p£ar  =  10%, 
i  =  1,  ...,11,  dL  <  d  <  du,1 
d  £  R14  and  X  £  R14 


Table  8  RBDO  Result  of  2-D  Mathematical  Example  with  Inaccurate  Design  Sensitivity 


Input 

Optimum  design 

Iter. 

Reli. 

Cost 

Probability  of  failure  (%) 

model 

d  1 

d  2 

cri 

02 

Anal. 

G\ 

G2 

G3 

A 

4.6310 

2.8821 

0.6946 

0.7205 

6 

20 

-1.3565 

3.5811 

2.2308 

0.5294 

B 

4.7572 

2.6993 

0.7136 

0.6748 

7 

18 

-1.4272 

2.6379 

2.2408 

0.5165 

Table  9  Performance  Measures  of  14-D  Engineering  Example 


Mode 

Function 

Value 

Feasibility  decision 

Full  frontal  impact 

G\ 

G2 

Chest  G 

Crush  displacement 

G3 

Brake  pedal 

G4 

Footrest 

Safety 

G5 

Left  toepan 

40%  offset  impact 

Ge 

Center  toepan 

<  Baseline  i 

Gi 

Right  toepan 

Gs 

Left  IP 

G9 

Right  IP 

NVH 

Gio 

Gu 

Torsion  mode 

Vertical  bending  mode 

1 0  Copyright  ©  20 1 5  by  ASME 

DISTRIBUTION  STATEMENT  A.  Approved  for  public  release;  distribution  is  unlimited. 


DISTRIBUTION  STATEMENT  A.  Approved  for  public  release;  distribution  is  unlimited. 


Table  10  Input  Model  for  14-D  Engineering  Example 


RVs 

Dist.  type 

Mean 

cov 

dB 

dL 

du 

Correlation 

Xi 

Normal 

d\ 

0.025 

1.9 

1.5 

2.3 

x2 

Normal 

d2 

0.025 

1.91 

1.5 

2.3 

X 3 

Normal 

di 

0.025 

2.51 

2.0 

3.0 

- 

x4 

Normal 

da, 

0.025 

2.4 

1.9 

2.9 

x5 

Normal 

ds 

0.025 

2.55 

2.0 

3.1 

x6 

Normal 

db 

0.025 

2.25 

1.8 

2.7 

Clayton  copula 

x7 

Normal 

dj 

0.025 

2.25 

1.8 

2.7 

r  =  0.5 

x8 

Gamma 

d  8 

0.02 

1.5 

1.2 

1.8 

Xio 

Gamma 

d<) 

0.02 

1.28 

0.9 

1.6 

X20 

Lognormal 

d  10 

0.025 

1.22 

0.9 

1.5 

X23 

Lognormal 

dn 

0.04 

0.75 

0.6 

1.0 

X25 

Lognormal 

dn 

0.04 

0.65 

0.5 

0.8 

Gaussian  copula 

X26 

Lognormal 

dn 

0.04 

0.85 

0.6 

1.1 

r  =  0.3 

XNi 

Weibull 

d\n 

0.04 

0.9 

0.7 

1.1 

- 

As  shown  in  Table  9,  there  are  11  constraints  in  this 
example.  They  are  safe  if  their  values  are  less  than  the  baseline 
values.  The  target  probability  of  failure  for  all  performance 
measures  is  10%.  Design  variable  vector  d  starts  from  the 
baseline  design  d5  and  has  to  be  in  between  lower  bound  AL  and 
upper  bound  dy.  Component  values  of  d5,  dL,  and  du  have  been 
shown  in  Table  10.  For  the  efficient  optimization  process, 
deterministic  design  optimization  (DDO)  is  performed  first. 
Generally,  the  DDO  optimum  design  is  close  to  the  RBDO 
optimum  design  because  the  cost  function  is  minimized  and 
active  constraints  are  determined  in  DDO.  The  RBDO  is  carried 
out  to  increase  the  reliability  (reduce  the  probability  of  failure) 
of  the  active  constraints  while  minimizing  the  cost.  Hence,  the 
RBDO  optimum  design  could  be  obtained  efficiently  by  starting 
RBDO  from  the  DDO  optimum  design.  The  sampling-based 
RBDO  has  been  carried  out  using  200,000  MCS  samples.  The 
optimization  process  is  completed  using  30  reliability  analyses. 
Among  the  30  analyses,  19  are  design  iterations,  which  means 
11  are  additional  line  searches.  This  could  be  a  good  indication 
that  the  provided  design  sensitivity  is  accurate  and  effective 
because  the  accurate  design  sensitivity  of  probability  of  failure 
reduces  the  total  number  of  reliability  analyses.  If  an  inaccurate 
design  sensitivity  was  used,  the  optimization  would  not  find  a 
better  solution  in  the  direction  of  the  design  sensitivity.  In  that 
case,  the  search  for  the  next  design  point  along  the  line  search 
direction  would  come  back  to  the  current  design  point,  and  all 
the  line  searches  would  be  wasted  in  that  design  iteration. 
Therefore,  the  entire  optimization  process  becomes  more 
efficient  when  an  accurate  design  sensitivity  is  provided.  As 
explained  earlier,  the  design  sensitivity  calculation  can  be  done 
during  the  reliability  analysis  (probability  of  failure  estimation) 


with  a  little  additional  cost.  Hence,  the  design  sensitivity 
calculation  itself  is  efficient  as  well.  Therefore,  the  developed 
design  sensitivity  method  would  provide  the  utmost  efficiency 
due  to  the  reduction  of  line  searches  and  its  efficient  calculation. 


Table  11  Optimal  Designs 


Design 

variable 

Initial 

design 

DDO 

RBDO 

d  1 

1.9 

1.8768 

1.8273 

d2 

1.91 

1.8857 

2.1891 

di 

2.51 

2.3008 

2.8540 

d\ 

2.4 

1.9753 

1.9701 

db 

2.55 

2.8003 

2.7060 

db 

2.25 

2.4358 

2.2611 

di 

2.25 

2.1066 

2.3265 

d  8 

1.5 

1.5353 

1.7841 

d<) 

1.28 

0.9 

0.9 

d  10 

1.22 

0.9 

0.9 

dn 

0.75 

0.6 

0.6 

d  12 

0.65 

0.5522 

0.5647 

d  13 

0.85 

1.1 

1.099 

d  14 

0.9 

0.7 

0.7 

The  obtained  optimum  designs  are  shown  in  Table  11.  In 
the  DDO  procedure,  d$~  dn  and  du  went  to  their  lower  bounds, 
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while  d  13  reached  its  upper  bound;  they  remain  on  their  bounds 
until  RBDO  converges.  The  other  variables,  d\  ~  dg  and  dn, 
have  changed  to  find  designs  that  are  more  reliable  than  the  DDO 
optimum  design  in  the  RBDO  procedure.  The  cost  value  and  the 
probabilities  of  failure  at  the  optimal  designs  are  shown  in  Table 
12.  In  the  DDO  process,  the  cost  function  is  greatly  reduced 
compared  to  the  initial  design  (baseline  design),  while  the 
probabilities  of  failure  are  not  reduced  as  much.  The 
probabilities  of  failure  are  reduced  below  the  10%  target 
probability  of  failure  after  RBDO;  at  the  same  time,  the  cost 
function  has  a  small  increment.  Again,  the  successful 
convergence  verifies  that  the  RBDO  is  performed  effectively 
using  the  developed  sensitivity  method.  Therefore,  it  can  be 
concluded  that  the  developed  design  sensitivity  method  performs 
well  for  RBDO  with  fixed  COV. 


Table  12  Cost  and  Probability  of  Failure  at  Optimal  Designs 


Design 

variable 

Initial 

design 

DDO 

RBDO 

Cost 

269.47 

257.58 

259.88 

G\ 

48.8% 

49.4% 

9.9% 

G2 

50.8% 

50.9% 

10.0% 

g3 

53.9% 

25.7% 

0.0% 

g4 

54.6% 

16.9% 

0.1% 

g5 

58.4% 

58.5% 

1.8% 

(j6 

58.8% 

65.4% 

9.9% 

Gi 

58.9% 

67.2% 

9.9% 

Gs 

53.3% 

55.9% 

9.9% 

G9 

51.5% 

52.3% 

6.9% 

G10 

49.3% 

0.1% 

0.0% 

Gn 

51.0% 

50.6% 

10.0% 

5.  CONCLUSION 

In  many  industrial  problems,  the  tolerance  of  an  input 
random  variable  is  specified  as  a  percentage  value  of  the  mean 
of  the  input  random  variable.  To  obtain  reliable  yet  less  costly 
design  in  the  problems,  RBDO  with  fixed  COV  should  be 
considered.  When  the  COV  of  an  input  random  variable  is  fixed, 
the  STDEV  changes  along  with  the  design  variable,  which  is  the 
mean  of  the  input  random  variable.  Then,  input  randomness 
fluctuates  due  to  the  change  of  the  STDEV.  As  a  consequence, 
the  optimum  design  of  RBDO  becomes  a  moving  target  problem 
due  to  the  fluctuation  when  COV  is  fixed.  Therefore,  accurate 
design  sensitivity  of  probability  of  failure  is  necessary  to  achieve 
a  correct  optimum  design  in  RBDO  with  fixed  COV. 

In  this  paper,  the  design  sensitivity  of  the  probability  of 
failure  has  been  developed  using  a  first-order  score  function  for 
sampling-based  RBDO  with  fixed  COV.  The  score  function 
contains  the  effects  of  changes  in  both  mean  and  STDEV. 


Therefore,  the  developed  design  sensitivity  can  consider  the 
change  of  the  STDEV  due  to  the  change  of  design  variable, 
which  is  the  mean  of  the  input  random  variable.  Correlated  input 
random  variables  are  considered  in  the  developed  design 
sensitivity  having  a  copula  density  function  in  the  score  function 
derivation. 

Using  a  highly  nonlinear  function,  the  developed  design 
sensitivity  of  the  probability  of  failure  has  been  tested  using  the 
FDM.  The  accuracy  of  the  developed  design  sensitivity  is 
verified  as  the  sensitivity  agrees  with  the  FDM  design  sensitivity. 
Moreover,  the  developed  method  uses  only  0.25%  MCS  samples 
of  the  FDM  design  sensitivity,  and  thus  quite  efficient.  The 
RBDO  with  fixed  COV  is  performed  for  a  2-D  mathematical 
example  using  the  developed  design  sensitivity.  The 
optimization  result  shows  smooth  convergence,  so  the 
effectiveness  of  the  developed  method  has  been  confirmed.  In 
addition,  for  the  same  example,  it  is  also  shown  that  ignoring  the 
effect  of  the  STDEV  in  the  design  sensitivity  causes  difficulty  in 
the  convergence  of  the  optimization.  Flence,  it  can  be  seen  that 
the  developed  design  sensitivity  is  necessary  for  RBDO  with 
fixed  COV.  Finally,  the  developed  method  is  applied  to  a  14-D 
engineering  example.  The  example  converges  in  19  design 
iterations  with  1 1  additional  line  searches  from  DDO  optimum 
design.  This  indicates  that  the  developed  sensitivity  provided 
very  accurate  information  in  the  optimization  direction  search 
and  that  the  accurate  design  sensitivity  helped  the  optimization 
process  use  few  design  iterations  so  the  whole  optimization 
process  could  be  efficient.  Flence,  it  is  confirmed  that  the 
developed  design  sensitivity  method  is  effective  and  efficient. 
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APPENDIX 


Table  A.1  Probability  Density  Function  (PDF)  and  Cumulative  Distribution  Function  (CDF)  of  Marginal  Distribution 


Distribution  type 

PDF  fx(x;  a,b) 

CDF  Fx(x\  a,b ) 

nx 

Normal 

1  1  (x-a\2 

— - e  2  \—J 

yJ2nb 

/X  —  CL\ 

*(  „  )* 

(— oo,  oo) 

Lognormal 

'I  l/ln  x-a\2 

r  ?\  h  1 

/In  x  —  a\ 

(0,  oo) 

y/2nxb 

b  ) 

Weibull 

b-(*S-'e-©b 

a  V a J 

(0,  oo) 

Gumbel 

1  r  X  —  Cl  /  X  —  CKn 

exp[  exp  (  b  )] 

(— oo,  oo) 

bexp[  b  exp  b  )\ 

Gamma 

g-x/b 

ra-l  a 

T(a)ba 

P(b,x/aT 

(0,  oo) 

Extreme 

1  r  X  —  Ci  (X  — 

rxpi  b  expU  )J 

1  exp  [  exp  (  b  )] 

(— oo,  oo) 

Extreme  Type  II 

“(b-)a+1e-& 
b  \xJ 

e-(l)a 

(0,  oo) 

Table  A.  2  Relationship  between  Mean,  STDEV,  and  Parameters  of  Marginal  Distribution 


Distribution  type 

Parameters 

Domain  of  parameters 

Normal 

li  —  a,  a  —  b 

a  E  (— oo,  oo),  b  E  (0,  oo) 

Lognormal 

H  =  ea+b2/2 

a2  =  (efa2  -  l)e2a+b2 

a  £  (— oo,  oo),  b  £  (0,  oo) 

Weibull 

(U  =  ar(l  +  1  /b) 

a 2  =  a2[r(l  +  2/b)  -  T2(l  +  1/A)] 

a  E  (0,  oo),  b  £  (0,  oo) 

Gumbel 

(tt  =  a  +  0.5772 b,  a2  —  b2n2/6 

a  £  (— OO,  oo),  b  £  (0,  oo) 

Gamma 

(U  =  ab,  a2  =  ab 2 

a  E  (0,  oo),  b  £  (0,  oo) 

Extreme 

li  —  a  —  0.5772 b,  a2  —  b2n2/6 

a  £  (— OO,  oo),  b  £  (0,  oo) 

Extreme  Type  II 

li  =  b  T(1  —  1/a) 

a 2  =  A2[r(l  —  2/a)  —  f2(l  —  1/a)] 

a  E  (0,  oo),  b  E  (0,  oo) 

*  $(•):  CDF  of  standard  normal  distribution 
§  H»:  Gamma  function,  T(s)  =  Jo°°  ts~1e~t  dt 

**  P(»,»):  Reguralized  gamma  function,  P(s,x )  =  y(s,x)/r(s)  =  J *  ts_1e_t  dt/P(s) 
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